Effects of HBT correlations on flow measurements 
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The methods currently used to measure collective flow in nucleus-nucleus collisions assume that 
the only azimuthal correlations between particles are those arising from their correlation with the 
reaction plane. However, quantum HBT correlations also produce short range azimuthal correlations 
between identical particles. This creates apparent azimuthal anisotropies of a few percent when 
pions are used to estimate the direction of the reaction plane. These should not be misinterpreted 
as originating from collective flow. In particular, we show that the peculiar behaviour of the directed 
^ ■ and elliptic flow of pions observed by NA49 at low pr can be entirely understood in terms of HBT 

' correlations. Such correlations also produce apparent higher Fourier harmonics (of order n > 3) of 

the azimuthal distribution, with magnitudes of the order of 1%, which should be looked for in the 
data. 
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I. INTRODUCTION 



In a heavy ion collision, the azimuthal distribution of particles with respect to the direction of impact (reaction plane) 
is not isotropic for non-central collisions. This phenomenon, referred to as collective flow, was first observed fifteen 
years ago at Bevalac jj], and more recently at the higher AGS and SPS || energies. Azimuthal anisotropies are 
very sensitive to nuclear matter properties ||J^]. It is therefore important to measure them accurately. Throughout 
■ this paper, we use the word "flow" in the restricted meaning of "azimuthal correlation between the directions of 
t-H | outgoing particles and the reaction plane" . We do not consider radial flow Q , which is usually measured for central 
0^ ■ collisions only. 

Flow measurements are done in three steps (see |7j for a recent review of the methods): first, one estimates the 
direction of the reaction plane event by event from the directions of the outgoing particles; then, one measures the 
azimuthal distribution of particles with respect to this estimated reaction plane; finally, one corrects this distribution 
for the statistical error in the reaction plane determination. In performing this analysis, one usually assumes that 
the only azimuthal correlations between particles result from their correlations with the reaction plane, i.e. from flow 
This implicit assumption is made, in particular, in the "subevent" method proposed by Danielewicz and Odyniec Q 
r* ■ in order to estimate the error in the reaction plane determination. This method is now used by most, if not all, heavy 
ion experiments. 

However, other sources of azimuthal correlations are known, which do not depend on the orientation of the reaction 
plane. For instance, there are quantum correlations between identicalparticles, due to the ( ant i) symmetry of the 
wave function : this is the so-called Hanbury-Brown and Twiss effect [[)[, hereafter denoted by HBT (see fnjjllj for 
reviews). Azimuthal correlations due to the HBT effect have been studied recently in In the present paper, 

we show that if the standard flow analysis is performed, these correlations produce a spurious flow. This effect is 
important when pions are used to estimated the reaction plane, which is often the case at ultrarelativistic energies, 
in particular for the NA49 experiment at CERN jl3). We show that when these correlations are properly subtracted, 
the flow observables are considerably modified at low transverse momentum. 

In section 2, we recall how the Fourier coefficients of the azimuthal distribution with respect to the reaction plane 
are extracted from the two-particle correlation function in the standard flow analysis. Then, in section 3, we apply 
this procedure to the measured two-particle HBT correlations, and calculate the spurious flow arising from these 
correlations. Finally, in section 4, we explain how to subtract HBT correlations in the flow analysis, and perform 
this subtraction on the NA49 data, using the HBT correlations measured by the same experiment. Conclusions are 
presented in section 5. 
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II. STANDARD FLOW ANALYSIS 



In nucleus-nucleus collisions, the determination of the reaction plane event by event allows in principle to measure 
the distribution of particles not only in transverse momentum px and rapidity y 7 but also in azimuth 0, where <f> is 
the azimuthal angle with respect to the reaction plane. The 4> distribution is conveniently characterized by its Fourier 
coefficients Q 

J cosn- A ' •' 
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dp 

where the brackets denote an average value over many events. Since the system is symmetric with respect to the 
reaction plane for spherical nuclei, (sinn</>) vanishes. Most of the time, because of limited statistics, v n is averaged 
over pt and/or y. The average value of v n (px, y) over a domain T> of the (j>t, y) plane, corresponding to a detector, 
will be denoted by v n (T>). In practice, the published data are limited to the n — 1 (directed flow) and n = 2 (elliptic 
flow) coefficients. However, higher harmonics could reveal more detailed features of the 4> distribution 0. 

Since the orientation of the reaction plane is not known a priori, v n must be extracted from the azimuthal correlations 
between the produced particles. We introduce the two-particle distribution, which is generally written as 

dN dN dN . 

(1 + C(pi,p 2 )) (2) 
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where C(pi,p 2 ) is the two-particle connected correlation function, which vanishes for independent particles. The 
Fourier coefficients of the relative azimuthal distribution are given by 

jOTcosnfai - fc) V^a Widfa 

Cn(PTl,yi,PT2,V2) = (COS n(0i - 2 )) = ~ JlT r 1 " • (3) 

a pict P2 

We denote the average value of c n over (j>T2, 2/2) in the domain T> by c u (pti, 2/1, "D ), and the average over both (pri, Vi) 
and {pT2,y2) by c n (V,V). 

Using the decomposition (||), one can write c„ as the sum of two terms: 

Cn(PTl,yi,PT2,y2) = cf™ (p T l , 2/1 , PT2 , V%) + c ^ on_flow (p T1 > y± ) p T2 , J/ 2 ) (4) 

where the first term is due to flow: 

Cn° W (pTl,yi,PT2,y2) = V „ (p T l , Vl)V n (pT2 , 2/2) (5) 

and the remaining term comes from two-particle correlations: 

JOT cos - 2 )C(pi , p 2 ) i N M d^dfo 
C n - fl ° W bTi,2/i,PT2,2/2) = — -™ d ~^P2 (6 ) 

tidwp-/^ 2 

In writing Eq.(|^), we have used the fact that (sinn</>i) = (sinn</> 2 ) = and neglected the correlation C(pi, p 2 ) in the 
denominator. 

In the standard flow analysis, non-flow correlations are neglected J^,^|, with a few exceptions: the correlations due 
to momentum conservation are taken into account at intermediate energies fl5|| , and correlations between photons 
originating from ir° decays were considered in |l6| ] . The effect of non-flow correlations on flow observables is considered 
from a general point of view in ]l7| . In the remainder of this section, we assume that c° on_flow = 0. Then, v n can 
be calculated simply as a function of the measured correlation c n using Eq.(|^), as we now show. Note, however, 
that Eq.(||) is invariant under a global change of sign: v n (pT,y) — * —v n (pr,y). Hence the sign of v n cannot be 
determined from c n . It is fixed either by physical considerations or by an independent measurement. For instance, 
NA49 chooses the minus sign for the v\ of charged pions, in order to make the v\ of protons at forward rapidities 
come out positive []l3| . Averaging Eq.(^|) over (pti,2/i) an d {PT2,V2) in the domain D, one obtains : 

v n (V) = ±y/c n (V,V). (7) 
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This equation shows in particular that the average two-particle correlation c n (X>,X>) due to flow is positive. Finally, 
integrating (JsJ) over (pt2 : IJ2), and using (Q), one obtains the expression of v n as a function of c„: 

v n (pT 1>yi ) = ± -■ (8) 

This formula serves as a basis for the standard flow analysis. 

Note that the actual experimental procedure is usually different: one first estimates, for a given Fourier harmonic to, 
the azimuth of the reaction plane (modulo 27t/to) by summing over many particles. Then one studies the correlation 
of another particle (in order to remove autocorrelations) with respect to the estimated reaction plane. One can then 
measure the coefficient v n with respect to this reaction plane if n is a multiple of to. In this paper, we consider only 
the case n = m. Both procedures give the same result, since they start from the same assumption (the only azimuthal 
correlations are from flow). This equivalence was first pointed out in [jl8|. 



III. AZIMUTHAL CORRELATIONS DUE TO THE HBT EFFECT 

The HBT effect yields two-particle correlations, i.e. a non-zero C(pi,p2) in Eq.(||). According to Eq.(|J), this 
gives rise to an azimuthal correlation c5J 011_fl ™, which contributes to the total, measured correlation c„ in Eq.([|). In 
particular, there will be a correlation between randomly chosen subevents when one particle of a HBT pair goes into 
each subevent. The contribution of HBT correlations to c5J°°~ flow will be denoted by c™ T . 

In the following, we shall consider only pions. Since they are bosons, their correlation is positive, i.e. of the same 
sign as the correlation due to flow. Therefore, if one applies the standard flow analysis to HBT correlations alone, i.e. 
if one replaces c n by c™ T in Eq.(||), they yield a spurious flow v™ T , which we calculate in this section. 

First, let us estimate its order of magnitude. The HBT effect gives a correlation of order unity between two identical 
pions with momenta pi and P2 if IP2 — Pi I ^5 h/R, where R is a typical HBT radius, corresponding to the size of 
the interaction region. From now on, we take h = 1. In practice, R ~ 4 fm for a semi-peripheral Pb-Pb collision 
at 158 GeV per nucleon, so that 1/R ~ 50 MeV/c is much smaller than the average transverse momentum, which is 
close to 400 MeV/c: the HBT effect correlates only pairs with low relative momenta. 

In particular, the azimuthal correlation due to the HBT effect is short-ranged : it is significant only if 2 — 4>i ^ 
\/{Rjpx) ~ 0.1. This localization in <f> implies a derealization in n of the Fourier coefficients, which are expected to 
be roughly constant up to n < Rpt ~ 10, as will be confirmed below. 

For small n and (pri ,yi) in 2?, the order of magnitude of c™ T (pti , Vi , "D ) is the fraction of particles in T> whose 
momentum lies in a circle of radius 1/R centered at pi. This fraction is of order (R 3 (pr) 2 '(ttit) Ay) -1 , where (px) 
and (tot) are typical magnitudes of the transverse momentum and transverse mass (tot — \/p\ + to 2 , where to is 
the mass of the particle), respectively, while Ay is the rapidity interval covered by the detector. Using Eq.(^), this 
gives a spurious flow of order 

I \ V2 



I"" BT(p) H «w( TO) aJ ■ < 9 > 

The effect is therefore larger for the lightest particles, i.e. for pions. Taking R = 4 fm, (px) ~ (tot) ~ 400 MeV/c 
and Ay = 2, one obtains |?; n (X>)| ~ 3 %, which is of the same order of magnitude as the flow values measured at SPS. 
It is therefore a priori important to take HBT correlations into account in the flow analysis. 

We shall now turn to a more quantitative estimate of c™ T . For this purpose, we use the standard gaussian 
parametrization of the correlation function (||) between two identical pions [l9[| : 

C Y (Pi,P 2 ) = \e-^ R >- q o R o-^ R2 L (10) 

One chooses a frame boosted along the collision axis in such a way that p\ z + p2z — ("longitudinal comoving 
system", denoted by LCMS). In this frame, qx, q and q s denote the projections of P2 — Pi along the collision axis, 
the direction of pi + P2 and the third direction, respectively. The corresponding radii Rl, R a and i? s , as well as the 
parameter A (0 < A < 1), depend on pi + p 2 . We neglect this dependence in the following calculation. Note that 
the parametrization (^) is valid for central collisions, for which the pion source is azimuthally symmetric. Therefore 
the azimuthal correlations studied in this section have nothing to do with flow. Note also that we neglect Coulomb 
correlations, which should be taken into account in a more careful study. We hope that repulsive Coulomb correlations 
between like-sign pairs will be compensated, at least partially, by attractive correlations between opposite sign pairs. 
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Since C(pi,p2) vanishes unless P2 is very close to pi, we may replace dN/d 3 p2 by dN/d 3 pi in the numerator of 
Eq.(^|), and then integrate over P2. As we have already said, q s , q and qL are the components of P2 — pi in the 
LCMS, and one can equivalently integrate over q s , q a and qL- In this frame, y\ ~ and one may also replace dN/d 3 pi 
by (l/rriTi)dN / d 2 pf ^yi. The resulting formula is boost invariant and can also be used in the laboratory frame. 

The relative angle <p2 ~ 4>i can be expressed as a function of q s and q Q . If pri 3> 1/i?, then to a good approximation 

02 - 4>i q s /pTi- (11) 
If pti ~ l/R, Eq. (|TT|) is no longer valid. We assume that R s ~ R and use, instead of ([ll]), the following relation : 

<?.? + 9o = Pti + Pt2 ~ 2 Pt\PT2 cos(</> 2 - <f>i). (12) 

To calculate c™ T (pn , yi , T>) , we insert Eqs.(|i"o|) and ( |TT| ) in the numerator of @ and integrate over (q s ,q ,qL). 
The limits on g D and <7l are deduced from the limits on (pt2,2/2), using the following relations, valid if pxi l/R : 



q a = PT2 - Pti 

q L = m T1 (y 2 - yi). (13) 



Since q s is independent of px2 and 2/2 (see Eq.(|ll|)), the integral over q s extends from —00 to +00. 

Note that values of q a and qL much larger than l/R do not contribute to the correlation ([lO]), so that one can 
extend the integrals over q a and qL to ±00 as soon as the point (pti, yi) lies in T> and is not too close to the boundary 
of D. By too close, we mean within an interval 1/R ~ 50 MeV/c in pr or l/(i?tmT) ~ 0.3 in y. One then obtains 
after integration 

1 dN 



hbt 1 iw _ Xn3/2 ( n2 \ m Ti d 2 pTidy 



(jPTi,yuV) = exp - ^ . (14) 

Jv d Pt 2 «2/2 



At low pt, Eq.(pTj) must be replaced by Eq.(p_2|). Then, one must do the following substitution in Eq.(|l 

where \ = RsPt and Ik is the modified Bessel function of order k. 

Let us discuss our result ([14]). First, the correlation depends on n only through the exponential factor, which 
suppresses c™ T in the very low px region pn ^5 n/2R s . For n smaller than R s (pr) — 10, the correlation depends 
weakly on n, as discussed above. Neglecting this n dependence, ( |l4| ) reproduces the order of magnitude (^|). To see 
this, we normalize the particle distribution in D in order to get rid of the denominator in (|l4|), and the numerator 
(l/mTi)(dN/d 2 pTidyi) is of order 1 / (pt) 2 (wt) Ay . However, Eq.(|T4|) is more detailed, and shows in particular that 
the dependence of the correlation on pxi and y\ follows that of the momentum distribution in the LCMS (neglecting 
the wit and y dependence of HBT radii). This is because the correlation c™ T is proportional to the number of 
particles surrounding pi in phase space. 

Let us now present numerical estimates for a Pb-Pb collision at SPS. We assume for simplicity that the pr and y 
dependence of the particle distribution factorize, thereby neglecting the observed variation of (pt) with rapidity [ p0[ . 
The rapidity dependence of charged pions can be parametrized by [po[ : 

' eJJl^m (16, 



dy cjVZk \ 2a2 , 
with er = 1.4 and (y) = 2.9. The normalized px distribution is parametrized by 

dN e m / T 



d 2 p T 27rT(m + T) 



exp 



with T ~ 190 MeV IgOJ. This parametrization underestimates the number of low-pr pions. The values of R , R s and 
Rl used in our computations, taking into account that the collisions are semi-peripheral, are respectively 4 fm, 4 fm 
and 5 fm The correlation strength A is approximately 0.4 for pions jp3[ . 

Finally, we must define the domain T> in Eq.([l4]). It is natural to choose different rapidity windows for odd and 
even harmonics, because odd harmonics have opposite signs in the target and projectile rapidity region, by symmetry, 
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and vanish at mid-rapidity ((y) — 2.9), while even harmonics are symmetric around mid-rapidity. Following the NA49 
collaboration Q, we take 4 < y < 6 and 0.05 < p T < 0.6 GeV/c for odd n, and 3.5 < y < 5 and 0.05 < p T < 2 GeV/c 
for even n. We assume that the particles in V are 85% pions half 7r + , half 7r~. Then, for an identified charged 
pion (a 7r + , say) with pt — pri and y — yi, the right-hand side of Eq.([l4|) must be multiplied by 0.85 x 0.5, which is 
the probability that a particle in T> be also a tt + . 
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FIG. 1. Apparent |un BT | from HBT correlations as a function of pr in MeV/c. Top : n = 1 (solid curve), n — 3 (long 
dashes), n = 5 (short dashes). Bottom : n = 2 (solid curve), n = 4 (long dashes), n — 6 (short dashes). The values are smaller 
for even n than for odd n because the rapidity intervals chosen to estimate the reaction plane differ. 

Substituting the correlation calculated from Eq.(|l4|) in Eq.(||), one obtains the value of the spurious flow v™ t (pt, y) 
due to the HBT effect. Fig.|l] displays |«™ T |, integrated between 4 < y < 5 (as are the NA49 data) as a function 
of pt- As expected, v^ BT depends on the order n only at low px, where it vanishes due to the exponential factor 
in Eq.(|l4|). HBT correlations, which follow the momentum distribution, also vanish if px is much larger than the 
average transverse momentum. Assuming that 1/R S <C rn,T, we find from Eq.(pd|) that the correlation is maximum 
at p T =PTmax where 



PTi 



IT 



T 



1/4 



nm 
2R S 



60v^ MeV/c, 



(18) 



which reproduces approximately the maxima in Fig.[j]. 

Although data on higher order harmonics are still unpublished, they were shown at the Quark Matter '99 conference 
by the NA45 Collaboration [^4| which reports values of v 3 and V4 of the same order as v\ and V2 , respectively, suggesting 
that most of the effect is due to HBT correlations. Similar results were found with NA49 data lE5fl. 



IV. SUBTRACTION OF HBT CORRELATIONS 



we can subtract this term from 
the measured correlation (left-hand side of Eq. (||) , which will be denoted by c m casurcd i n this section) to isolate the 
correlation due to flow. Then, the flow v n can be calculated using Eq.(^J), replacing in this equation c„ by the corrected 
correlation c„ ow = c™ oasured — c™ T . In this section, we show the result of this modification on the directed and elliptic 



flow data published by NA49 for pions 13 



The published data do not give directly the two-particle correlation c measured^ rather the measured flow u™ easured . 
Since these analyses assume that the correlation factorizes according to Eq.(|^), we can reconstruct the measured 
correlation as a function of the measured v n . In particular, 

measured/ <t-\\ measured/ \ measured /<t->\ /in\ 

c n \PTl,Vl,V) — V n (PTl,yi)V n (V). (19) 

We then perform the subtraction of HBT correlations in both the numerator and the denominator of Eq.(^). 

The integrated flow values measured by NA49 are v ^^(p) = -3.0 ± 0.1% and u^ured^) = 3.0 ±0.1% |lj. 
After subtraction of HBT correlations, the coefficients are smaller by some 20% : vi(T>) — —2.5% and V2(T>) — 2.6%. 
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Fig.g displays the rapidity dependence of v\ and v-i at low transverse momentum, where the effect of HBT correla- 
tions is largest. Let us first comment on the uncorrected data. We note that w ™ easured [ s ze ro below y < 4 (i.e. outside 
T>, where there are no HBT correlations) and jumps to a roughly constant value when y > 4 (where HBT correlations 
set in). This gap disappears once HBT correlations are subtracted, and the resulting values of v\ are considerably 
smaller. The values of Vi are also much smaller after correction, except near mid-rapidity. 




FIG. 3. Directed flow v\ and elliptic flow vi of pions, integrated between 4 < y < 5, as a function of the transverse momentum 
Pt in MeV/c, measured by NA49 (u™ casurcd , short dashes) and after subtraction of HBT correlations (v n , full curves). The 
long dashes show linear and quadratic fits at low pr for t>i and V2 , respectively. 

Fig.|| displays the pr dependence of v% and i>2- The behaviour of v n {pT) is constrained at low px : if the momentum 
distribution is regular at px = 0, then v n (pT) must vanish like pVp. One naturally expects this decrease to occur on a 
scale of the order of the average pr- This is what is observed for protons ]l3| ] . However, the uncorrected l; m easurcd an d 
^measured £ Qr pj ong rema j n large far below 400 MeV/c. In order to explain this behaviour, one would need to invoke a 
specific phenomenon occurring at low px- No such phenomenon is known. Even though resonance (mostly A) decays 
are known to populate the low-pr pion spectrum, they are not expected to produce any spectacular increase in the 
flow. 

HBT correlations provide this \ow-px scale, since they are important down to 1/R ~ 50 MeV/c. Once they are 
subtracted, the peculiar behaviour of the pion flow at low pt disappears. v\ and vi are now compatible with a 
variation of the type v\ oc pt and vi cx p\, up to 400 MeV/c. 



V. CONCLUSIONS 



We have shown that the HBT effect produces correlations which can be misinterpreted as flow when pions are 
used to estimate the reaction plane. This effect is present only for pions, in the (pt,u) window used to estimate the 
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reaction plane. Azimuthal correlations due to the HBT effect depend on px and y like the momentum distribution in 
the LCMS, i.e. (l/mx)dN/dyd 2 px, and depend weakly on the order of the harmonic n. 

The pion flow observed by NA49 has peculiar features at low px'- the rapidity dependence of v\ is irregular, and 
both v\ and V2 remain large down to values of px much smaller than the average transverse momentum, while they 
should decrease with pt as px and p T , respectively. All these features disappear once HBT correlations are properly 
taken into account. Furthermore, we predict that HBT correlations should also produce spurious higher harmonics 
of the pion azimuthal distribution (v n with n > 3) at low px, weakly decreasing with n, with an average value of 
the order of 1%. The data on these higher harmonics should be published. This would provide a confirmation of the 
role played by HBT correlations. More generally, our study shows that although non-flow azimuthal correlations are 
neglected in most analyses, they may be significant. 
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